Negative effective mass transition and anomalous transport in power-law hopping 

bands 
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We study the stability of spinless Fermions with power law hopping Hij oc |i — jl'" ■ It is shown 
that at precisely ac = 2, the dispersive inflection point coalesces with the band minimum and the 
charge carriers exhibit a transition into negative effective mass regime, m* < characterized by 
retarded transport in the presence of an electric field. Moreover, bands with a < 2 must be accompa- 
nied by counter-carriers with m* > 0, having a positive band curvature, thus stabilizing the system 
in order to maintain equilibrium conditions and a proper electrical response. We further examine 
the semi-classical transport and response properties, finding an infrared divergent conductivity for 
1/r hopping (a = 1). The analysis is generalized to regular lattices in dimensions d = 1, 2, and 3. 

PACS numbers: 71.10.Ca, 72.10.-d, 71.55.-1, 75.30.Hx 

I. INTRODUCTION 

Conducting models with long-ranged, power law tight- 
binding bands (PLTB) have been introduced in a number 
of different physical scenarios such as in DNA chains^, 
Frenkel excitons-, and resonating valence bond (RVB) 
phases'^. Recent theoretical work on PLTB bands has 
focused on disorder-induced critical phenomena, where 
critical delocalization occurs at a = 1 for random 
hopping^i^ and a mobility edge in the regime 1 < a < 3/2 
for non-random hopping and diagonal disorder^. Ad- 
ditionally, others have investigated variational magnetic 
states in strongly correlated Hubbard models at large 
onsite repulsion U >> t resulting in the t — J approxi- 



First let us clarify as to which energy regimes these 
band curvature effects play a significant role. The two 
primary limits to consider are near fc — > and k = kp. 
The former case is near the bottom of the band where the 
concavity and the functional form of the dispersion re- 
lation will govern the semi-classical dynamics and trans- 
port in the presence of electric and magnetic fields. How- 
ever, in the latter case these curvature effects are irrele- 
vant near the Fermi surface, where any dispersion can be 
expanded as follows (in ID for example). 
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mation with both long-ranged exchange J^, 



and hoppingiiiSiii t, 
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Earlier numerical inves- 



tigations of non-disordered PLTB systems by Borland 
and Mencherc® explored non-extensive effects of PLTB 
bands. They found that for finite chains of length N , 
for < a < 1 , there is anomalous wave-packet spreading 
and for a > 3/2, ballistic motion is recovered as the 
nearest-neighbor limit is approached. However, there 
have been no systematic studies of the dispersion rela- 
tions, or the transport properties of the non-disordered 
free-field model. 

In this article we study various properties of PLTB 
chains as a function of the parameter a. Here we re- 
strict our analysis to those particular cases of a that 
are tractable analytically, resulting in simple closed form 
representations of the dispersion relations. The effective 
mass is examined as a function of a and it shown pre- 
cisely that at ac = 2, the system enters into a negative 
effective mass regime, characterized by retarded trans- 
port in the presence of an electric field. Mathematically, 
this occurs at a critical value ac where the dispersive 
inflection point moves to the bottom of the band until 
the long wavelength behavior is dominated by a negative 
curvature. For uncoupled isotropic tight-binding bands, 
our results and methods can be generalized to regular 
lattices in dimensions 1,2, and 3. 



and the linear contribution is usually retained as e(fc) ~ 
vp^k — kp) near the Fermi points ±kp without any loss 
of generality. Therefore, irrespective of the curvature 
away from ep, the low-energy Fermi-surface instabilites 
should remain unaffected as long as the dimensionality 
and Fermi-surface topology remains unchanged. 



II. THE NON-INTERACTING BAND 
STRUCTURE 

A. The dispersion relations 

Consider the tight-binding model in d=l (the analysis 
to higher dimensional isotropic lattices) with one orbital 
per site, for N sites with lattice spacing a, governed by 
the following Hamiltonian (neglecting spin). 



(2.1) 



where ipj and ipj are fermion fields that obey ^f/'lj V'il 



The effective width of the spectrum tij has the power 



law spatial dependence. 
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where the constant depends on the details of the 
atomic orbital overlap matrix elements. For the entirety 
of this article we shall focus on the specific even positive 
integers a = 1,2 and 4, given that the full dispersion 
relations can be treated analytically, rather than only 
retaining asymptotic forms near the top and bottom of 
the band. Moreover, this particular range of exponents 
encompasses the earlier ranges of observed criticality in 
disordered systems. Because as \i — j\ oo, a 

proper convergence of the sums taken for a large chain is 
ensured. For translationally invariant systems with peri- 
odic boundary conditions we can make use of the Bloch 
wave representation. 
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where the fc's run over the first Brillouin zone. In the 
limit of a large system size the band dispersion takes the 
form, 
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FIG. 1: A comparison of the kinetic dispersion relations at 
different hopping powers of a {h — l).(a) At a = 1, there is a 
logarithmic divergence at A; = 0. (b) For the case a — 2 cor- 
responding to inverse-squared hopping, a Dirac cone emerges 
at smaller values of k unlike the Tomonaga-Luttinger liquid, 
which is linearized around the two Fermi points ikp, rather 
than at fc — * 0. (c) At shorter ranged hopping, the conven- 
tional parabolic band is restored, with some additional cur- 
vature, (d) Numerical determination of the inflection point 
I{a) as a function of a, where at precisely Oc = 2, /(oc) = 
(summations taken for N=500). 



where Lia{z) is the polylogarithm function^. For the 
aforementioned cases of a, one can make use of the follow- 
ing exactly summable series, which are valid for \k\ < 27r, 
which is beyond the natural cutoff of the 1st Brillouin 
zone^(Note that we have absorbed a factor of 2 into to): 
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The Eas. (|2.6ll2.8p have been plotted in Fig[T] Apparently, 
the preceding expressions produce the full band curvature 
within the first Brillouin zone, however in certain cases 
the explicit periodicity is absent but can be repeated to 
construct higher Brillouin zones. Time reversal invari- 
ance is preserved such that Ek.a — £-k,a- Note that the 
analysis can generalized to all isotropic regular lattices 
by simply reproducing each dispersion independently in 
each direction. For example the a = 1 case in cZ = 3 
cubic lattice can be written as. 
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relations, near fc — > 0. The first and second derivatives of 
ek,a respectively give rise to the group velocity Va{k) ~ 



1 c)efc.o 



and the effective mass ra*^{k) = ^- 



Taking the appropriate derivatives we have for the group 
velocity. 
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and the effective mass. 
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which have been plotted in FigO It follows that an in- 
fiection point I{a) can be defined as the location where 
m^ik) = 0. 



III. THE NEGATIVE EFFECTIVE MASS 
REGIME 

A. Semi-classical dynamics 



Apparently, the range of the interactions has a notice- 
able effect on the precise functional form of the dispersion 



In the usual particle/hole theory, electrons with charge 
e < are accelerated from the bottom of a parabolic 
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FIG. 2: The group velocity and effective mass of the electrons 
{h — 1) (a) Q = 1, the velocity diverges as the effective mass 
tends to zero for small k. (b) At a = 2 the velocity steadily 
decreases for a constant negative effective mass, (c) a = 4 
is relatively short-ranged, with a velocity that increases then 
decreases, with a shift to negative effective mass values orig- 
inating from an inflection point in the dispersion, (d) The 
effective mass for each case of a. Notice the emergence of 
negative values at a = 1,2 for small k, where the inflection 
point has disappeared, and the curvature of the dispersions 
becomes purely negative. 



B. Transition at q = 2 

For a large range of a's considered, there is a definite 
transition from a positive to negative band curvature. By 
only considering the exact expression given by Eas. (|2.1ip 
it is clear that 2 < ac < 4. Take the extreme case 
of short ranged hopping, which is the nearest-neighbor 
limit and can be approached at large a, and the dis- 
persion becomes e„„(A:) = —2^0 cos A; with m* > near 
the band minimum, having an inflection point exactly 
at /„„ = ±7r/2. The other limit of interest, can be ex- 
amined more precisely by performing the summations of 
Eq. (|2.5p numerically as shown in FiglT] (d) . Notice that 
as a becomes smaller, the inflection point moves closer to 
the band minimum, I{ac) = until the curvature has re- 
versed sign, at a « 2. Because the summations converge 
rapidly the finite size effects are negligible. Addition- 
ally, for a more precise analysis, one can make use of the 
derivative relation ^^^gjj^ = Lis_i(e^). Consequently, 
the effective mass is controlled by the function Lis{z) at 
negative values of s when a < 2 and in order to discern 
the behavior of the Lis{z) function at negative indices, 
we can employ the expansion^, 

^ m 

ii-.W = 7Y^^w+TE«™^fc^' (s>0) (3.2) 



band when an external positive electric field is applied. 
Electrons are therefore defined to have a positive effective 
mass, TO* > 0. On the other hand, holes possess the 
opposite charge e > and move in a (valence) band, with 
a negative curvature, generating a negative effective mass 
TO* < and thus will move in the opposite direction as 
the electron. Hence, the concavity and curvature of a 
band, will govern the sign and magnitude of to* > as a 
function of k. Semi-classically, the Boltzmann equation 
states that the fields E and H will change the local carrier 
concentration fk at the ratei^, 
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Consequently, electrons and holes, with reversed disper- 
sive concavity do not have the conventional metallic re- 
sponse to an applied electric field; there is no acceleration 
of the carriers, rather the opposite occurs and the elec- 
tronic transport is suppressed by a current /<^^. If such 
carriers participate in the conduction process, then a pos- 
itive current recorded by the measuring instrument must 
be accompanied by an internal current lo = + ly, 
where J> is part of the positive current due to particles 
having m* > 0. Thus, multiple band with both types 
of carriers should be introduced, such that a net positive 
effective mass can maintain equilibrium. 



where the coefficients are Eulerian numbers that satisfy 
the recurrence relations, 

am,k = (to + 1 - k)am~iM~i + kam-i,k (3.3) 

Taking the appropriate form of Eq. (|2.5p , it is clear that 
TO* < in the range < a < 2. 

C. Transport in 1/r bands 

For the longest ranged hopping exponent considered 
(a = 1), there is a unique logarithmic divergence at long 
wavelength values of fc, which is accompanied by a van- 
ishing effective mass, as shown in Figs. [IJa) and[2Ka), 
and with the reduced inertia the velocity is divergent. 
Such a singularity will have important consequences and 
requires a careful interpretation if such an equilibrium 
state is valid. The semi-classical conductivity tensor in 
the relaxation-time approximation is given by—, 
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where the integrations are taken over the Fermi surface, 
T{ep) is the scattering lifetime at the Fermi level and 
M is the effective mass tensor. For systems with cubic 
symmetry at arbitrary filling, the tensor reduces to a 
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scalar yielding 
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which is logarithmically divergent at the infrared limit, 
which coincides with the previously observed anomalous 
wavepacket spreading^, now also exhibiting anomalous 
transport characteristics. Another point to consider is 
the behavior of the Fermi level Ep as a function of the 
filling fraction or density. At a = 1 the Fermi energy 
depends on the density p as. 
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Upon inspection, Ea. (j3.6p vanishes at p = 1/3 or a one- 
third filling fraction unlike the nearest-neighbor band 
which has a Fermi energy that vanishes at half filling 
OT p = 1/2. In general charge/spin density wave forma- 
tion, superconductivity and Mott metal-insulator transi- 
tions are strongly filling fraction dependent. Therefore, 
it would be of interest to determine if there are p = 1/3 
band instabilities in the system in the presence of inter- 
actions, as a negatively diverging Fermi energy indicates 
that the a — 1 system is unstable at lower densities. 



D. 1/r^ hopping 

Of much interest is close to the transition point of 
a = 2, or inverse squared hopping. It follows that at 
small fc, the dispersion is linear, £fc_2 oc fc, which maps 
onto a pseudo-relativistic kinetic energy, analogous to 
the 6 Dirac points of graphene. However, unlike the 
Tomonaga-Luttinger liquid, ek,2 is not linearized close to 
the two Fermi points ±/cf and does not require a cutoff 
in order to remedy negative energy states^. It should be 
noted that other models that are relevant to a = 2 hop- 
ping were studied by Haldane^ and Shastrj*^, who deter- 
mined the ground state of the spin 1/2 antiferromagnetic 
Heisenberg chain with 1/r^ exchange. The Haldane- 
Shastry model was exactly solved by a Gutzwiller pro- 
jected wavefunction, identifying this particular class of 
inverse squared exchange models with Anderson's RVB 
phaseiS,. For a PLTB band the exchange constant J, in 
the limit of large repulsion U » t, is given by 
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Thus, the point here to emphasize is that long range 
exchange and hopping are related, as it requires a 1/r 
kinetic band to produce a 1/r^ exchange constant and 
therefore the earlier 1/r and negative effective mass 
transport characteristics become obviously important. 
Later, Kuramoto et al^ have included inverse squared 
hopping in addition to inverse-squared exchange, in a 



supersymmetric t — J model, showing that it is also 
exactly solved by a Gutzwiller projected wavefunction. 
Hence, the a ~ 1,2 cases with negative effective mass 
carriers are relevant to particular classes of RVB phases 
and strongly correlated systems. The consequences of 
this apparent connection needs to be explored further as 
the long-ranged hopping appears to control the magnetic 
frustration. 



E. Bulk properties for q < 2, d=l 

Let us further investigate the properties of the neg- 
ative effective mass regime in ID. Another method for 
probing the stability of the ground-state is to study the 
total energy. The emergence of a complex or imaginary 
component of the thermodynamic ground-state energy 
suggests that the system is in a non-equilibrium state. 
The density for a degenerate Fermi system is p = kp/i: 
and the total kinetic energy per particle can be obtained 
by integrating over the internal energy distribution. 
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The density of states can by easily computed by invert- 
ing the dispersion relations in terms of k and using the 
relation g{X) ~ 2/nhv{X). After substituting the velocity 
expressions from (j2.10p . we have 
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Subsequently, deep in the negative effective mass regime, 
for a — 1, the total energy becomes, 
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where is a complex constant. Ui is purely real as 
a function of density in the range (0, 1), indicating that 
there are no obvious bulk instabilities. Next let us exam- 
ine the ground-state energy at a = 2, which is in the crit- 
ical regime. After substituting Eg. p.lOp into Eq. (|3.8p . 
we have the following density p dependent ground-state 
energy 
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which is purely real, indicating a stable ground-state 
along with a stable Fermi energy. Lastly, the low energy 
scaling in terms of the density is often used to compare 
with Coulomb interactions or disorder. The well-known 
competition between the kinetic (band) energy and other 
terms can result in a change of the ground-state. There- 
fore, a naive low energy scaling of the total energies re- 
sults in the following, 



Uiip) cc Hp) 

U2{p) OC p 



(3.13) 
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where Ea. (|4.4p exhibits interference fringes(or the math- 
ematical form of beating) instead the usual case of 
purely oscillatory behavior shown in Eq. (|A7|) . Moreover, 
the condition for constructive interference requires that 



Q 



2kp)x = 27rn for n = 0,1,2.... This suggests 

that in inverse-squared hopping conducting materials, 
the umklapp momenta can be derived from the inter- 
ference fringe maxima of the Friedel oscillations. We can 
further examine the special case of the half-filled band 

where Q — Akp- This leads to 



IV. FERMI SURFACE NESTING AND 
SCREENING AT a = 2 



Evidently, carriers with m* < 0, will screen a charged 
impurity differently from the free Fermi gas. Therefore 
it would be useful to examine the susceptibility of the 
particle-hole channel to highlight some of these differ- 
ences. Lindhard theory is an important perturbative ap- 
proach to quantify the change in the electronic charge 
density dn due to a static impurity potential 0^^^*. In 
Appendix [X] we present the conventional formalism used 
for the ID free Fermi gas. Here, we shall focus on the 
analytically tractable ID, T=0 behavior of the static sus- 
ceptibility Xo{<l) at the negative mass transition(a — 2). 
Subsituting the form for ek,2 into Eq. ()A2|) yields the fol- 
lowing. 
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(4.1) 

which unlike the parabolic band or nearest neighbor case, 
does not exhibit a nesting singularity exactly at q — 2kF, 
rather it is shifted to q — ±2{tt — /cj?). It follows that we 
can restore the lattice constant and it becomes clear that 
the nesting is lattice dependent, as the spacing of the 
reciprocal lattice determines the nesting. 
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This dependence has futhcr implications for the RKKY 
function, which for a delta function impurity potential is 
simply the Fourier transform Eq. (|4.ip . 
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Consequently the at large distances we have, 

S2n(x) oc — - cos ({ Q — 2kF)x\ — cos f ( Q 
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which contains unusual Akp periodicity, reminiscent of 
the ID Wigner crystal density-density oscillations^. 



V. CONCLUSION 

In summary we have explored various properties and 
the general dependence of the band curvature on the hop- 
ping exponent a for power-law banded electron systems. 
Our methods and approach can be generalized to alterna- 
tive lattices in higher dimensions, and this article should 
serve as a basis to do so. We have demonstrated, using 
a variety of approaches, that the curvature completely 
changes sign in the range ac = 2 and have highlighted 
various static properties of the a = 2 system, including 
unusual lattice dependent Fermi surface nesting prop- 
erties. In addition, our investigations are particularly 
useful for those studying Anderson localization and criti- 
cal phenomena in Hamiltonians that contain long-ranged 
band structures as previous studies on such systems have 
highlighted range 1 < a < 1.5 as a critical regime. Here 
have demonstrated that even in the absence of disorder, 
the a = 1 system is not stable at lower densities and 
contains logarithmic divergence of its Fermi energy and 
anomalous transport characteristics. 

In closing we emphasize that the curvature of a band 
is generally associated with the magnitude and sign of 
the effective mass, resulting in important implications 
for electronic devices that require non-equilibrium condi- 
tions. For example, it has been shown that for negatively 
charged carriers with a negative effective mass, the cur- 
rent density is negative and the both the Hall emf and the 
thermal emf are positive, which follows directly from the 
Lorentz force law^^. Therefore, the capability of tuning a 
conducting system into a negative effective mass regime 
could be exploited for technological purposes, given that 
the interplay of non-equilibrium conditions with disorder 
and quantum/thermal fluctuations can affect the crit- 
icality and the transport properties. Additionally, an 
unexplored connection with strongly correlated models 
utilizing PLTB bands leading to magnetically frustrated 
spin-liquid states will be explored in a future work. 
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APPENDIX A: STATIC RESPONSE OF A ID 
FREE FERMI GAS 

The general dielectric response e{q) can be evaluated 
from the following linear response relation, 



i-0^"*(«)x(9) 



(Al) 



where 4>{q) is the full physical potential in q space and 
the static Lindhard response function in ID is given by, 
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where V denotes the principal part of the integral and fk 
is the Fermi-Dirac distribution. For the conventional ID 
non-interacting electron gas at T = 0, this evaluates to 
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which has the well known singularity at q = 2kp, associ- 
ated with Fermi surface nesting and spin/charge density 



wave instabilities. Furthermore, regarding the screening 
behavior, the simplest case assumes a delta function im- 
purity potential of the form. 
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which leads to the following real space density modula- 
tion: 



(A5) 



After subsitituting Eq. (jA3P , one has the following expres- 
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which is the well known Ruderman-Kittel-Kasuya-Yosida 
(RKKY) range function. The large distance asymptotic 
becomes, 



cos(2A:fx) 



(A7) 



displaying the 2k p periodicity and 1/x decay commonly 
referred to as Friedel oscillationsii. 
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